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Abstract 

We study a modified non-linear Schrodinger equation on a 2 dimensional sphere with radius R 
aiming to describe electron-phonon interactions on fullerenes and fullerides. These electron-phonon 
interactions are known to be important for the explanation of the high transition temperature of 
superconducting fullerides. Like in the i? — > cxd limit, we are able to construct non-spinning as well 
as spinning solutions which are characterised by the number of nodes of the wave function. These 
solutions are closely related to the spherical harmonic functions. For small R, we discover specific 
branches of the solutions. Some of the branches survive in the — > oo limit and the solutions 
obtained on the plane {R = oo) are recovered. 
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I. INTRODUCTION 



Carbon exists in several forms in nature. One is the so-called fullerene which was dis- 
covered for the first time in 1985 FuUerenes are carbon-cage molecules in which a big 
number n of carbon (C) atoms are bonded in a nearly spherically symmetric configuration. 
These C„ configurations typically have diameters of(i~7 — IsA and consist of 12 pentagons 
and — 10) hexagons, where > 24 has to be even. The most prominent example is the 
Cgo fullerene which also is called a "buckminster" fullerene. It consists of 20 hexagons and 
12 pentagons and has a diameter of d^o ^ 7 A. In these fuUerenes, three of the valence 
electrons of the C-atom are used to form the bounds with the neighbouring C-atoms, while 
the remaining "free" valence electron can "hop" along the n positions in the fullerene. 

Due to unoccupied electron orbitals, it is easy to reduce fuUerenes, i.e. put extra electrons 
on them. It was found that if alkali metal atoms-which donate one electron each- such as 
rubidium (Rb), potassium (K), caesium (Cs) or sodium (Na) are put onto a Cgo fullerene, 
this leads to a metallic or even superconducting behaviour [2]. The transition temperature 
Tc of these alkali-doped fuUerenes, also called fullerides, is rather high for superconductors, 
e.g. ^ 33 K fa a RbC..C,„ g The .upe.conduct.vity can be explained by phonon- 
electron interactions in the fullerides jjj. In [5| it was found that if a Cgo or a C70 molecule 
is doped with one or two excess electrons, the additional charges accumulate nearly along 
an equatorial line of the molecule. The distortion in the lattice and the electron distribution 
was found to be polaron-like. 

Recently, a modified non-linear Schrodinger equation has been studied in the context 
of fuUerene-related structures, so-called nanotubes 0, Q. These objects are graphite-like 
sheets roled to form a cylinder. The carbon atoms are bonded in a hexagonal lattice on this 
sheet. The existence of solitons which are created through the interaction of the "excitation" 
described by a complex scalar field ip and the lattice vibrations was discussed. The creation 
of these type of solitons was introduced by Davydov in the 1970s to explain the dispersion 
free energy transport in biopolymers. ip then is an amide I-vibration. Here, however, we 
think of ip as the wave function of an "excess" electron interacting with the lattice via 
electron-phonon interactions. The electron could e.g. be thought of as the valence electron 
of an alkali metal atom with which the fullerene is doped. The full discrete system of 
equations (which results from a Frohlich type Hamiltonian) can be approximated by the 
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above mentioned modified non-linear Schrodinger equation on the 2-dimensional plane. 

In this paper, we study the modified non-linear Schrodinger equation for the hexagonal 
lattice on a 2-dimensional sphere, thus aiming to describe the electron-phonon interactions 
in fullerenes and fullerides. These interactions are, as mentioned previously, of particular 
interest for the superconductivity of the fullerides. Our model is, of course, an approxima- 
tion since fullerenes consist of hexagons and pentagons. However, since a) there are more 
hexagons in the Ceo and b) the coefficients appearing in the non-linear Schrodinger equa- 
tion don't depend strongly on the actual form of the lattice, we beheve that this is a good 
approximation. 

Finally, let us mention that studying solitons on a sphere can appear to be useful also 
from a mathematical point of view, as was stressed e.g. in Q]. In the present case this 
appears to be true as well. In particular, we notice a strong relationship between most of 
the various branches of solutions and the spherical harmonics. 

Our paper is organised as follows: in Section II, we present the modified non-linear 
Schrodinger equation on the sphere, in Section III, we discuss non-spinning solutions and in 
Section IV spinning generalisations. We give our conclusions in Section V. 



II. MODIFIED NON-LINEAR SCHRODINGER EQUATION ON THE SPHERE 

We consider the following modified non-linear Schrodinger equation in 2 dimensions 
[lo| . which in dimensionless variables reads: 

i^ + A^ + agip (iV-r + feAl^p) = (1) 



where a, h are constants with a = 2, h = ^ for the square lattice studied in 10 1 and 

n 

a = 4, 6 = g for the hexagonal lattice studied in '6]. The additional term proportional to 
h results from both the discreteness of the lattice and the interaction between the lattice 
and the complex scalar field %p ( "phonon-electron" interaction), g is the ( "phonon-electron" ) 
coupling constant of the system and determines the "strength" of the non-linear character 
of the equation. Especially, as can be seen from the derivation of the modified non-linear 
Schrodinger equation from the discrete equations Q|, the mass of the carbon atom is encoded 
into g. 

As stated in the introduction, we will study the above equation on a 2-dimensional sphere 
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with radius R. The Laplacian operator then reads: 

1 / (9^ d I d"^ \ 

^ = ^(^5^ + -t^5^ + -^^J • (2) 

As usual, the angles are such that e [0 : tt] and e [0 : 27r] and ip = tp{t, 6, if) is a complex 
valued scalar field. Since we are mainly concerned with fullerenes here which consist (next 
to pentagons) out of hexagons, we use in analogy to Q| a = 4 and & = |- The solutions of 
(0) can be characterised by their norm rj: 

r^^ = R^ j \tp\^ sin edOdip (3) 

as well as by their energy: 

E = R^ [ - -g\^\^ + 2bg{V\ij\Y) smOdOd^ . (4) 



2- 

The coupling constant g can be absorbed into the scalar field ip by rescaling it as follows: 
ip '4'/y/9- III our numerical simulations, we will set g = 1. This in turn means that 
we drop the normalisation constraint [rj = 1) for the wave function ip. However, since we 
present explicitely the norm of our numerical solutions, the normalized wave function and 
the corresponding value of g can be easily computed. Especially, if a solution ceases to exist 
for a specific value of the norm, this means that the soliton solution exists only above a high 
enough value of the phonon-electron coupling. 

III. NON-SPINNING SOLUTIONS 



To construct explicit solutions, we use the following axially symmetric Ansatz ll|, 3| ■ 



^ = e''^'^{9) . (5) 

Inserting the Ansatz into equation (^, we find: 

+ cot e<^' + + ^f^i?^ = (6) 

where the prime denotes the derivative with respect to 6. Note that we are fixing the length 
scale A of the space variable by Rphy = XR, where Rphy is the physical radius of the fuUerene 
and A is of the order of A. We will vary R here in order to be able to take the R oo limit 
..s ccnpa. ou. w.h those of gQ. TKe otK. di— fu. pKysica. ,„anti«es 



can be recovered by reinserting the physical constants of the model appropriately (see e.g. 
eq. (6), (7) of ^j). As far as the frequency Uphy is concerned, we have: 

where m is the mass of the "electron". Thus, since u is dimensionless, ujphy has dimension 
of an inverse time: [uJphy] =sec~^. 

The regularity of the solutions at 5 = 0, tt requires the following boundary conditions : 

$'|e=o = $'|e=. = (8) 

For later convenience, we denote by A4 the reflexion operator with respect to the equatorial 
(9 = n/2) plane : M^{9,ip) = $(7r - 9, if). Solutions for which = $ we will call 

symmetric, those with = — $ antisymmetric and finally solutions for which ^ ±$ 
will be called asymmetric in the following. 

It appears impossible to construct explicit solutions of (jHI), (©• This is why numerical 
techniques have to be adopted. However, in our numerical study, it became obvious that the 
pattern of solutions occuring for different values of u is quite complicated. In an attempt to 
understand this pattern, we have analysed the linearized equation first. This analysis will 
be discussed prior to presenting the numerical results. 

A. Linearized equation 

First, we remark that (jB)) has three fixed points, repectively $c = and = ^s/^- 
Perturbating around one of these fixed points according to 

^9) = ^, + efi9) + Oie') (9) 

and linearizing the equation in e leads to 

f" + cot9f -ujR^f = if $^ = (10) 

and 

In both cases, we recognize the equation of Legendre polynomials P„ : P" + cot 9P^ + 
n(n + l)Pn = 0. If we take into account the boundary conditions given by (jH)), the relevant 



solutions of the linearized equation are given in terms of these polynomials and fix a relation 
between R, uj and the integer n: 

f{e) = Pn{e) with ujR^ = -n{n + 1) if = (12) 



and 



f{e) = PM with 



2uR^ 



n{n + 1) if $c = ± 




(13) 



In particular, the above relations define two sets of critical values of the spectral parameter 
as a function of R which are indexed by n. For (jl2|) and (jl3|) we find respectively : 



From the above analysis, setting n = 1 such that Pi oc cos 6', we conclude that the 
pertubated solution about $c = a/co'/4 is positive (i.e. especially it has no nodes, at least 
for values of uj close to the critical values) and has a maximum above the north pole (6' = 0), 
a minimum above the south pole {9 = tx). 

The perturbated solution about $c = leads to a solution which is antisymmetric under 
the reflexion M.. It vanishes at the equator (i.e. has a node at 9 = 7r/2) and possesses a 
maximum above the north-pole, a minimum above the south pole. 

Of course, the functions presented here constitute only approximations (even up to an 
unfixed multiplicative constant) to the solutions of the full equation (jHl). The closer we are 
to the critical points (fT^ - (fTH|l . the more accurate these approximations are. Nevertheless, 
preforming this analyis we gain extremely useful information about the solutions of the full 
equations, namely about (i) their symmetries and (ii) the critical values of u where the 
various branches start or finish. In the following, we will discuss the numerical results. 

B. Numerical results 

We constructed solutions numerically for generic values of R and of uj. In our numerical 
routine, $(0) is used as a shooting parameter and the regular solution as well as the cor- 
responding value of UJ are determined numerically. We first constructed a "fundamental" 
solution which has no node {k = 0). For generic values of the parameters, this solution 



n{n + 1) 



(14) 



i?2 




(15) 
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corresponds to the above mentioned deformation of the Legendre polynomial Pi and is thus 
positive, has a maximum for ^ = and a minimum at = tt. This is illustrated in Fig.l 
for $(0) = 1.0, 0.5 and 0.3. Keeping R fixed and varying $(0) we obtain a branch of asym- 
metric solutions labelled by uj. Our numerical analysis strongly suggests that these regular 
solutions exist for uJcr,i < uj < oo with 

- • (16) 



In the limit tu CJcr,i the solution tends to the constant solution $ = A/ct'cr,i/4, which for 
i? = 2 is equal to 0.2582. As can be seen in Fig.l, the solution for $(0) = 0.3 is already 
nearly constant and thus close to this critical solution. 

In correspondence with these solutions, the mirror symmetric ones also exist and consti- 
tute another branch with the same values of u. 

Next, we constructed a family of solutions having one node for 6 G]0,7r[. The node 
occurs at 6 = 7i/2 and the solutions of this branch are antisymmetric under the reflection 
Ai, as illustrated in Fig. 2. These solutions exist only for uj > uJcr,i = and for 

arbitrary large values of uj. In the limit uj —>■ uJcr,i, the solutions are of the form with 
$c = and f{6) = cos6. This is shown in Fig. 2, where the solution with $(0) = 0.25 
is proportional to the trigonometric function cos^. When uj is increased the effects of the 
non-linear term become stronger and the solution's profile deviates considerably from an 
elementary trigonometric function. 

When investigating this branch of solutions for larger values of $(0) (or equivalently of uj) 
we observed that the one-node solution develops a plateau (surrounding the region 9 = 7r/2) 
on which the function $ is very close to zero. Away from this plateau the solution resembles 
two disconnected solitons located at the two poles of the sphere. This is illustrated in Fig. 
3 (solid line) for $(0) = 3.7 and the corresponding value of uj is uj = 15.95. However, 
for large enough values of $(0), i.e. of uj, our numerical analysis strongly indicates that 
asymmetric solutions exist next to the antisymmetric ones described above. In Fig. 3 such a 
asymmetric solution corresponding to uj = 14.85 together with a comparable antisymmetric 
solution is presented. The asymmetric solution in Fif. 3 (dotted line) has an energy of 
order -900 (in our dimensionless variables). This has to be contrasted with the energy of the 
antisymmetric solution (solid line) which is of order -300. This result is interesting with view 



to 



B. 



There, it was found that the excess electron density in the Cqq fuUerene is located 



along an equatorial line of the Ceo- The excess electron density here can be represented 
by $^(6'). We find here, that the solution with the electron density located at the equator 
(asymmetric solution) is energetically more favourable than the one which has the electron 
density located in one of the hemispheres. It is remarkable that the two solutions show a 
very similar behaviour around 6 = and start to deviate from each other at roughly 6 > 0.8. 
Unfortunately, this similarity of the solutions in the region around the north pole renders the 
numerical construction of the solutions extremely difficult. This is why we did not succeed 
to construct the full asymmetric branch. However, let us mention that a similar feature 
occurs as well in the spinning case = 1. In this case, the numerical analysis was much 
easier and we will thus discuss this point again in the section about spinning solutions. 

The values of the norm rj, the energy E and the value $(0), which corresponds to the 
maximal value of the solution, are shown in Fig. 4 as functions of uj for both the fundamental 
{k = 0) and the one-node {k = 1) solution with R = 2 (the graphic is limited to u < 1.5). 

Solutions having more nodes of the function $(^) can be constructed in correspondence 
with the higher order Legendre polynomials. However, the systematic study of these solu- 
tions is not the aim of this paper. 

Finally, we remark that 



This is in agreement with |7[ , and confirms that in the plane both fundamental and one- node 
solutions can be constructed for a; > 0. 

IV. SPINNING SOLUTIONS 

In order to construct spinning solitons, we use the Ansatz Q| : 



lim u!, 

R^oo 



'cr,l — 



lim UJ, 



'cr,l — 



. 



(17) 




(18) 



which, after inserting into (0), leads to the equation : 




1 + 5f<l>2 sin^ 9 



1 



. 



(19) 



Regularity of the solutions at 6* = 0, vr requires 



$(0) = $(7r) = . 



(20) 
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In this case, $c = is the only fixed point of the equation and the corresponding hnearized 
equation is solved by the associated Legendre functions , —n < N < n, provided the 
relation (|14p holds. In what follows, n will be related to the number k of nodes of the 
solution hj k = n — N. 

We would like to mention that we were able to construct an exact solution of the above 
equation for = 1 and i?^ = 3/4. It is of the form 

$ = a sm 6* with u = — , (21) 

3 

where a; is a real constant. Obviously, this solution is symmetric. 
A. Numerical results for A'^ = 1 

In this case, the numerical solutions can be constructed by using the value $'(0) as a 
shooting parameter. We studied numerically the solutions corresponding to the deformation 
of the associated Legendre function oc sin 6 (for no-node solutions) and of P2 oc sin 6 cos 9 
(for one-node solutions) and found that in accordance with (fT^ . there exist no solutions for 
respectively 

2 6 

UJ < = UJcr,l , UJ < -— = UJcr,2 (22) 

As for the case of non-spinning solution, the norm rj, the energy E and the parameter $'(0) 
are presented in Fig. 5. 

As can be seen from this figure, a new phenomenon occurs in the case of spinning solutions 
without nodes. Namely, the main branch corresponding to the deformation of stops 
at some critical value of uj, say ujmax with Umax > ^cr,i (for P = 2, we find uj^ax ~ 0.3). 
However, the inspection of Fig. 5 by no means reveals that some critical phenomenon occurs 
in this limit. When plotting the quantities P, uj and the norm rj as functions of the shooting 
parameter $'(0) (see Fig. 6), though, it is clearly indicated that some critical phenomenon 
appears. Indeed as our numerical results indicate, the slope of the curves becomes infinite 
when the shooting parameter reaches the critical point. The solutions on this branch are 
symmetric under the operator M.. The profile of the symmetric solution corresponding to 
$'(0) = 0.25 is shown in Fig. 7. 

For sufficiently large values of tu, say uj > uj, ( U ^ 1.0 for P = 2) a new branch of 
nodeless solutions appears which seems to exist for arbitrary large values of a;, while for 
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iJ > io > ujmax no solutions seem to exist at all. 

The solutions for uj > uJ are asymmetric under the reflexion Ai . Only the limiting solution 
a.t uj = uJ is symmetric under Ai. An asymmetric solution corresponding to $'(0) = 1.0 as 
well as the limiting, symmetric solution for $'(0) = 0.36 is shown in Fig. 7. Clearly the 
asymmetry increases with increasing $'(0). In other words, when u > uJ increases, the 
wave function of spinning, nodeless solitons is more and more concentrated in one of the 
hemispheres. 

All the above results were obtained for R = 2. However, we also studied the dependence 
of the critical values of u on the radius R of the sphere. The results are given in Table 1. 



Table 1: Critical values of uj for different R 



R 




^max 


UJ 


oo 











5 


-2/25 


0.03 


0.15 


2 


-1/2 


0.3 


1.0 


1.5 


-2/(1.5)2 


1.3 


2.15 


v/3/4 


-3/8 


oo 


oo 



For generic value of R such that a/3/4 < R < oo two branches of solutions exist, namely 
a branch of symmetric solutions for u G [uJcr,i '■ ^max] and a branch of asymmetric solutions 
for G [cJ : oo]. In the limit R oo , the symmetric branch disappears and only the 
asymmetric one survives. In the limit R a/3/4 only symmetric solutions seem to exist 
since u; ^ oo. In this respect, the explicit solution (PT|) plays a major role. 

Finally, we constructed one branch of spinning, one-node solutions. They exist for 
u! > — 6/i?2 = uJcr,2 and are symmetric under the operator Ai. They vanish when uj 
approaches the lower critical limit ujcr,2 and seem to exist for arbitrary large values of uj. 
Spinning, one node solutions corresponding to $'(0) = 1.0, 2.0 and 3.0 are shown in Fig. 8. 
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B. Numerical results for N = 2 



We have also constructed spinning solutions corresponding to = 2. The solutions are 
deformations of the associated Legendre functions P2 (xl — cos(20) (for no-node solutions), 
respectively P| oc cos 9 — cos(30) (for one- node solutions). In agreement with the analyis of 
the linearized equation, no solutions exist respectively for 

6 12 . , 

^ < = t^cr,2 , < = t^cr,3 • (23) 

Furthermore, the no- node solutions are symmetric under Ai, while the one-node solutions 
are antisymmetric under Ai. No branch of asymmetric solutions seems to exist in this case. 
Thus, as far as our numerical analysis indicates, in contrast to the = 1 case, these are 
the only N = 2 solutions for n < 3. 



V. CONCLUSION 

The study of the modified non-linear Schrodinger equation of |6|, |7|, llO[ on the sphere 
reveals a very rich set of stationary solutions. They can be characterised by their behaviour 
under the reflection operator Ai with respect to the equator {6 = tt/2) plane. We constructed 
branches of solutions which are either i) symmetric under A4, ii) antisymmetric under Ai or 
iii) asymmetric under A4. The solutions are further characterised by their spinning number 

and the number of nodes k of the wave function which we find -as expected- to be closely 
related to the quantum numbers of the spherical harmonics. Linearizing the modified non- 
linear Schrodinger equation, we find that the lower value of u, above which solutions exist, 
is determined by the spectral parameters obtained from the linearized equation. 

An unexpected phenomenon occurs in the cases N = 0, k = 1 and N = 1, k = 0. For 
A^ = 0, = 1, we find that for comparable values of uj two branches of solutions exist: one 
which is antisymmetric and one which is asymmetric. For A^ = 1, = 0, we observe that 
for finite R > a/3/4 the symmetric solutions only exist on a finite interval of the parameter 
UJ. Since uj is related to the kinetic energy of the solution [7], this means in other words 
that only for solutions with sufficiently low kinetic energy is the probability density of the 
"electron" the same in both hemispheres. If the kinetic energy becomes large enough, the 
electron density is essentially concentrated in one of the hemispheres. In p, it was found 
that the excess electron density in fuUerides is located around the equator. Here, we find 
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$^(^) to be concentrated around the equator in two cases: (a) for the N = 0, k = 1 case if 
u is large enough, (b) ior N = 1, k = ii u E [LJcr,i{R) '■ i^max(-R)]- 

We have not studied N > 2 and/or /c > 1 in this paper. We beheve that the above 
mentioned properties hold also true for these type of solutions and that, in particular, the 
solution with A^, k is related to the associated Legendre function P^,j>^- 
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FIG. 1: The profile of the non-spinning, nodeless {N = 0, k = 0) solution and its derivative is 
presented for i? = 2 and three values of ^*(0), namely $(0) = 1.0, 0.5 and 0.3. 
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FIG. 2: The profile of the non-spinning, one-node {N = 0, k = 1) solution and its derivative is 
presented for i? = 2 and three values of ^*(0), namely $(0) = 1.0, 0.5 and 0.25. 
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FIG. 3: The profile of one solution on the antisymmetric branch (corresponding to u = 15.95) as 
well as one solution on the asymmetric branch (corresponding to a; = 14.85) is shown for N = 0, 
k = 1 and R = 2. The corresponding derivatives of the solutions are also shown. 
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FIG. 4: The values of the energy E, the norm r] and of $(0) are shown as functions of u for the 
fundamental {k = 0) and one-node {k = 1) solutions with N = (non-spinning) and R = 2. 
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FIG. 5: The values of the energy E, the norm r] and of $'(0) are shown as functions of u for the 
fundamental {k = 0) and one-node {k = 1) solutions with N = 1 (spinning) and R = 2. 
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FIG. 6: The values of the energy E, the norm r/ and of cu are shown as function of $'(0) for the 
spinning, nodeless {N = 1, k = 0) solutions with R = 2. 
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FIG. 7: The profile of tlie spinning, nodeless {N = 1, k = 0) solution and its derivative is shown 
for = 2 and three values of $'(0), namely $'(0) = 1.0, 0.36 and 0.25. 
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FIG. 8: The profile of the spinning, one- node {N = 1, k = 1) solution and its derivative is shown 
for i? = 2 and three values of $'(0), namely $'(0) = 3.0, 2.0 and 1.0. 
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